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It is well-known that the Haldane phase of one-dimensional spin-1 chain is a symmetry protected 
topological (SPT) phase, which is described by a nonlinear Sigma model (NLSM) with a Q— term 
at = 2-7T. In this work we study a 3+1 dimensional SPT phase of SU(2JV) antiferromagnetic spin 
system with a self-conjugate representation on every site. The spin ordered Neel phase of this system 
has a ground state manifold M = ^jvjxutiv) ' ant ^ ^ ms system is described by a NLSM defined with 
manifold M. Since the homotopy group 7r4[.M] = Z for N > 1, this NLSM can naturally have a 
0— term. We will argue that when Q — 2ir this NLSM describes a SPT phase. We will also construct 
a 3d lattice model for this SPT phase, whose long wavelength physics is precisely described by this 
NLSM with = 2ir. This SPT phase is protected by the SU(2iV) spin symmetry, or its subgroup 
SU(iV) xSU(iV) xi Z2, without assuming any lattice symmetry. 

PACS numbers: 



1. INTRODUCTION 

In quantum many-body systems, quantum mechan- 
ics can lead to at least three types of exotic/nontrivial 
quantum disordered phases: (1) Topological phases with 
a gapped spectrum and bulk topological degeneracy, 
(2) algebraic liquid phases with gapless bulk spectrum 
and power-law correlations, and (3) symmetry protected 
topological phases. 

A symmetry protected topological (SPT) phase is a 
state of matter with gapped and nondegenerate bulk 
spectrum, but cannot continuously evolve into a direct 
product state without a bulk phase transition, when and 
only when the Hamiltonian of the entire evolution is in- 
variant under certain global symmetry G [l[ . In terms of 
its phenomena, a SPT phase on a d— dimensional lattice 
should satisfy at least the following three criteria: 

(1) . On a d— dimensional lattice without boundary, 
this phase is fully gapped, and nondegenerate; 

(2) . On a d— dimensional lattice with a (d — 
1)— dimensional boundary, if the Hamiltonian of the en- 
tire system (including both bulk and boundary Hamilto- 
nian) preserves certain symmetry G, this phase is either 
gapless, or gapped but degenerate. 

(3) . The boundary state of this d— dim system cannot 
be realized as a (d — l)-dim lattice system built with the 
same onsite Hilbert space, and with the same symmetry 
G. 

If a d— dim SPT phase satisfies all three criteria (1), (2) 
and (3), this phase is a SPT. Both the 2d quantum spin 
Hall insulator and 3d Topological insulator 

are perfect examples of SPT phases protected by time- 
reversal symmetry. 

Notice that the second criterion (2) implies the follow- 
ing two possibilities: On a lattice with a boundary, the 
system is either (2a) gapless, or (2b) gapped but degener- 
ate. When d > 3, the degeneracy of (26) can correspond 



to either spontaneous breaking of G, or correspond to cer- 
tain topological degeneracy at the boundary. Which case 
occurs in the system will depend on the detailed Hamil- 
tonian at the boundary of the system. For example, with 
interaction, the edge states of 2d QSH insulator, and 3d 
TBI can both be gapped out through spontaneous sym- 
metry breaking at the boundary, and this spontaneous 
symmetry breaking can occur through a boundary tran- 
sition, without destroying the bulk state 

In this work we will focus on bosonic spin systems. 
The simplest example of SPT phase of spin system is the 
Haldane phase of one dimensional spin-1 chain. Unlike 
the free fermion case, although there is a classification of 
bosonic SPT using group cohomology [l| , specific models 
of higher dimensional bosonic spin systems are not well 
understood. So far, in most studies, construction of 2d 
and 3d bosonic SPTs has been focused in systems with 
U(l) symmetry 11-13]. In this work we will study a 3 + 1 
dimensional analogue of the Haldane phase, which is con- 
structed as a SU(2iV) spin state with a self-conjugate rep- 
resentation on each site. Just like the Haldane phase, this 
3+ld SPT is described by a nonlinear sigma model de- 
fined with semiclassical antiferromagnetic order parame- 
ter plus a topological O— term. 



2. HALDANE PHASE 

Although symmetry protected topological phase is a 
pure quantum phenomenon without any classical ana- 
logue, the Haldane phase can still be described semiclas- 
sically by a nonlinear Sigma model (NLSM), which is 
defined only in terms of the semiclassical Neel order pa- 
rameter n [3, mi : 

C = - (d^n) 2 + ^-e abc e^n a d l ,n b d„n c . (1) 
a »7T 
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When the system has S0(3) symmetry, the entire mani- 
fold of the configurations of Neel order parameter is S 2 . 
If we assume the trivial vacuum has = 0, then when 
= 2tt, this model describes the Haldane phase. At the 
boundary of the Id system, i.e. the domain wall between 
O = and 2tt, the 0— term reduces to a O+ld Wess- 
Zumino-Witten (WZW) term at level 1, whose ground 
state has two fold degeneracy [lfj. As long as the SO(3) 
symmetry is preserved, this two fold degeneracy at the 
boundary cannot be lifted. Thus by definition the Hal- 
dane phase is a SPT protected by SO (3) spin rotation 
symmetry. 

Now let us couple two Haldane phases to each other: 

1 i2n 
£ = -(df.rii) 2 + ——tabc^^nld fl n\d v n c 1 
9 on 

+ 1^2 + A(Hi-n 2 ). (2) 

When A = — oo, effectively n\ = n2 = ft, then the system 
is effectively described by an 0(3) NLSM with O = 4n; 
while when A — +oo, the effective NLSM for the sys- 
tem has = 0. When parameter A is tuned from — oo 
to +oo, the entire phase diagram with A G (— oo,+oo) 
is gapped in the bulk. Thus the theory with O = 4n 
and = are equivalent. This analysis implies that 
with SO(3) symmetry, Id spin systems have two differ- 
ent classes: there is a trivial class with O = 4irk, and a 
nontrivial Haldane class with O = (4k + 2)ir. This 
classification is consistent with the cohomology formal- 
ism developed in Ref. [lj. I n fact, using the renormaliza- 
tion group calculation [17H20| . and the general nonper- 
turbative argument in Ref. |2l| , we can derive a phase di- 
agram for model Eq. [T] The system is topological when 
O G ((4k + 1)-7T, (4k + 3)7r), while the system is trivial 
when G ((4k - 1)tt, (4k + l)vr); © = (2k + I)tt is the 
transition, where the system is either gapless, or two fold 
degenerate. 

There are various ways of describing the Haldane phase 
on the lattice. In the follows we will choose one particular 
description that will be generalized to higher dimensions 
later. The Haldane phase can be described on the lattice 
as follows: On every site we introduce a slave fermion 
with both spin-1/2 index and SU(2) color index: /i,A,ai 
and the spin-1 operator is represented as [22j 

1 2 

^ ' = fj.A.a^apfj.A.^- (3) 

Z A=l 

In order to match the slave fermion Hilbert space with 
the spin-1 Hilbert space, we have to impose two different 
constraints on each site: 

^2 fi,A,afaA,a = 2 > ^2 f},A, a PAB-f 0,B ,a = 0, (4) 
a, A a,A,B 

where p M are three Pauli matrices of the color space. The 
second constraint guarantees that on every site the color 



space is in a total antisymmetric singlet, thus the spin 
is in a total symmetric spin-1 representation. Then the 
Haldane phase corresponds to the following mean field 
state of /a, a- fi,a forms valence bonds on links (2i, 2i + 
1), while f2, a forms valence bonds on links (2i+ 1, 2i + 2) 
(Fig. [TJi) . In terms of low energy field theory of the slave 
fermion, the Haldane phase is described by the following 
Lagrangian: 

C = ■07 P (9 M '0 + m i -ii)p z ip. (5) 

Here the Dirac fermion ip is the low energy mode of /, 
which is expanded at the two Fermi points kf = ±7r/2 
in the Id Brillouin zone. If we couple the Neel order 
parameter to the slave fermion, 

(-!)% ■ fj,A, a ^pfj,A,p ~ ft ■ lfa5Vll>, (6) 
A 

Eq. [T] can be readily derived after integrating out the 
slave fermions [23[ • Notice that in Eq. [5] the gauge fields 
introduced by constraints Eq. [4] are ignored, but in 1+1 
dimension gauge fields are always confining. 

The Haldane phase is a SPT only when the color- 
singlet constraint fjP^fj = is strictly imposed on every 
site, i.e. when the Hilbert space on every site is rigorously 
spin-1. If this constraint is given up, the Hilbert space 
on every site is enlarged to 6 dimension, and the Haldane 
phase becomes trivial, because it can now be adiabati- 
cally connected to a direct product state with spin-singlet 
on every site. Actually, besides the Haldane phase mass 
gap m4, we can consider another mass gap m5^75/9 2 V' 
of the Dirac fermion ip. Physically 7715 corresponds to a 
"color density wave" on the lattice, which is not allowed 
if the color singlet constraint is imposed strictly on every 
site. Without the color singlet constraint, the Haldane 
phase can adiabatically evolve into the color density wave 
state, by turning on 7715. 

Eq. Q] already assumed a rigorous spin-1 Hilbert space, 
since it does not involve any orbital degree of freedom 
besides vector n. 

3. 3D SPT OF SU(2A0 SPIN SYSTEM 
3.1 Field Theory Description 

Let us try to look for higher dimensional generaliza- 
tions of the Haldane phase of spin-1 chain, which has a 
description in terms of NLSM plus a O— term. The most 
naive generalization would be the AKLT state in higher 
dimensions, for instance the spin-2 AKLT phase on the 
square lattice. The boundary of the spin-2 AKLT phase 
on the square lattice is a spin-1/2 chain, which according 
to the LSM theorem cannot be gapped and nondegenrate, 
thus the spin-2 AKLT state seems to be a SPT. However, 
in order to protect the spin-2 AKLT state, we need trans- 
lation symmetry, since the boundary spin-1/2 chain can 
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FIG. 1: (a). The pictorial representation of the Haldane 
phase, where the blue and red links stand for the valence 
bonds of slave fermions with color A — 1 and 2 respectively. 
(6) A domain wall of 0' on the boundary of our 3d SPT. (c) 
The Young diagram of the self-conjugate representation of the 
SU(2iV) spin system that we are considering. 



be dimerized and gapped out once the translation sym- 
metry of the system is explicitly broken. Thus this is not 
an ideal generalization of the Id Haldane phase, whose 
stability does not rely on any translation symmetry. The 
spin-2 AKLT state on the square lattice can indeed be 
described by a NLSM with a topological term, but the 
configurational space of this NLSM would involve both 
the Neel and dimerization order parameters. 

The goal of this paper is to find a three dimensional 
SPT phase without assuming the translation symmetry. 
Inspired by Eq. [TJ we should first look for magnetic sys- 
tems, whose ground state manifold M. of the spin ordered 
phase has a nontrivial homotopy group 7T4 [M] — Z, then 
this SPT can be described by a NLSM defined in man- 
ifold M with a 9 term. The SU(2iV) antiferromagnet 
with self-conjugate representation satisfies this criterion: 
its magnetic ordered phase has GSM [H, [25| 



M = 



U(2AT) 



U(A0 x U(iV)' 



ir 4 [M] = Z, for N > 2. (7) 



Every Neel order configuration V G M can be repre- 
sented as 

liVxJV, Onxn \ 
V = VtaV, tl = I (8) 
Onxn, -Ijvxjv / 

where V is a SU(2iV) matrix. V is a hermitian traceless 
order parameter that satisfies V 2 — 1. In fact, when N — 
1, M. is precisely S 2 , and V can always be represented 
as V — n ■ a, where n is the Neel vector. 

Since n 4 [M] = Z, the following NLSM defined on M 
can be written down: 



S = J d 3 xdr haldpVdpV] 



256tt : 



■tr[Pd^Vd v Pd p VdxV}e^ pX . 



(9) 



By tuning the parameter g, there is obviously an order- 
disorder transition. When g is small, the system is in 
a spin ordered phase where V is condensed and sponta- 
neously breaks the SU(2A^) symmetry; when g is large, 
the system is in a disordered phase, and this disordered 
phase is what we are interested in. 

In the follows we will focus on the disordered phase 
of Eq. [S] with = 2ir, while assuming the trivial vac- 
uum of this spin system has = 0. With = 27r, the 
bulk spectrum of the field theory is identical to = 0, 
thus the disordered phase is gapped and nondegenerate. 
Then we can safely integrate out the bulk, and look at 
the boundary theory. Using the standard bulk-boundary 
correspondence, we can derive the boundary theory of 
Eq. |9j which is a 2+1 dimensional NLSM defined in M 
with a WZW term at level k = 1: 



Sb= J d 2 xdr ^[d^Vd^P] 



du I d 2 xdr 



256tt 2 



tT[Vd^Vd v Vd p VdxP]e ll J e \Q) 



Here u € (0, 1), and V[x, r, u) is an extension of V(x, r) 
that satisfies 



V{x,t,1)=V(x,t), V(x,T,0) = n. 



(11) 



Such WZW terms can be analyzed very reliably in 
0+ld and 1+ld, and in both cases, these terms change 
the ground state dramatically. In 0+ld, a WZW term 
leads to degenerate ground states; in 1+ld, it drives the 
system to a stable gap less fixed point described by con- 
formal field theory 2^, 27 1. In higher dimensions, non- 
trivial effects of a WZW term are still expected, but we 
no longer have a complete understanding. In the follows 
I will argue that the disordered phase of the 2+ Id bound- 
ary theory Eq. [TOlmust be either gapless or degenerate. 

In order to make this argument, let us first 
weakly break the SU(2iV) symmetry down to 

SU(A0xSU(A0xZ 2 . This residual SV(N)xSV(N) 
symmetry transformation can be written as 



V L , 



NxN 



V = 



OjVxTV, Vr 



(12) 



while the residual Z<i symmetry corresponds to exchang- 
ing Vl and Vr. With this symmetry reduction, the order 
parameter V can be written as 



V 



cos(0)1at x jv, ism(8)U 
— i sin(#)J7 1 ', — cos(9)ljy X N 



(13) 



where U is an SU(A r ) matrix. Under the transforma- 
tion Vl and Vr, U transforms U — > vIuVr, thus the 



4 



SU(iV) xSU(iV) residual symmetry precisely corresponds 
to the left and right transformation of U. Under the Z 2 
symmetry transformation, 



17 -> uK 



(14) 



With this symmetry reduction, the boundary theory 
with the WZW term Eq. [TU] is reduced to a principal 
chiral model (PCM) defined on manifold SU(7V) with a 
0' term: 



S b -+ J d 2 xdr hrld^d^U] 

2AlT 



(15) 



If the Z 2 symmetry of the residual symmetry is unbroken, 
i.e. the expectation value of 9 is tt/2, then the derived 
boundary SU(iV) PCM has precisely 9' = jr. When 
this Z 2 symmetry is explicitly broken, namely at the 
boundary we turn on a background field that tunes 9 
away from tt/2, then a straightforward calculation shows 
that the derived 0' is also tuned away from jr, and 
9' = 27r(2 + cos(60)(sin((9/2)) 4 . Under the Z 2 transfor- 
mation, 6' -> 2?r - 0'. 

The phase diagram of 2+ld SU(iV) PCM was studied 
in Ref. [2l|, where it was argued that when g is large 
enough, 0' > 7r and 0' < tt are two different disor- 
dered phases, which are both fully gapped and nonde- 
generate. These two disordered phases are separated by 
either a first or second order transition at 0' = jr. Thus 
at 0' = 7r the system must be either gapless or two fold 
degenerate. In our formalism we can see that the residual 
SU(-/V)xSU(-/V)xZ 2 symmetry guarantees 0' = 7r at the 
boundary, i.e. these symmetries guarantee the boundary 
of Eq. IH1 cannot be trivially gapped out. 

Just like the Id case discussed before, the existence of 
a ©'—term in a 2+lc? PCM can lead to nontrivial edge 
states. Since the 2d boundary of a 3c? system cannot 
have its own boundary, the edge states can only exist at 
a domain wall of physical order parameters. Now let us 
create the following configuration of 0' at the 2d bound- 
ary (Fig.QJ): 



&{x) = 2tt, for \x\ < R, 



&{x) = 0, for |x| > R. 



(16) 
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in the disordered 
R. there is a 



According to our previous work 
phase of Eq. [T5l at the domain wall \x 
gapless SXJ(N)l x S\J(N)r conformal field theory with 
level k = 1. The SU(iV) L and SU(N) R charges move 
clockwise and counter-clockwise respectively along the 
domain wall. Later these domain wall states will help 
to argue that with the SU(7V)xSU(A r )xZ 2 symmetry, 
model Eq. [T5] can only be realized at the boundary of a 
3d system, i.e. Eq. [9] describes a 3d SPT. 



By coupling two copies of Eq. [9] together like Eq. [2l 
one can show that the theory with = Att can be con- 
tinuously connected to = without a bulk transition. 
Thus again Eq. [9] describes a 3d SPT with Z 2 classi- 
fication: there is a trivial class with = Ajrk, and a 
nontrivial Haldane class with = (Ak + 2)tt. Using the 
same argument as Ref. 2l[ , we can derive a phase dia- 
gram for model Eq. |9j The system is topological when 
E ((Ak + l)jr, (Ak + 3)7r), while the system is trivial 
when 6 ((Ak - 1)tt, (Ak + l)?r); = (2k + l)jr is the 
transition, where the system is either gapless, or two fold 
degenerate. 



3.2 Lattice Construction 

Now let us construct a spin state for Eq. Hon a lat- 
tice. Consider a SU(2iV) antiferromagnet on a diamond 
lattice, where there are two flavors on each site, and the 
spin of each flavor on every site carries a self-conjugate 
representation of SU(27V) (Fig. QJ;). The spin operator 
can be represented using slave fermion fi t A,a,a : 



2 1 



A=l 



1 

2N 



Sal: 



(17) 



Here A — 1,2 is a color index, a, (3 — 1,---2N are 
the SU(2iV) indices, a = 1,2 denotes the two flavors. 
In order to match the spin Hilbcrt space and fermion 
Hilbcrt space, again we need to impose two constraints 

T,A, a fi,A, a ,Ji,A,a,* = 2JV, and fl A p%f l .B = 0, whose 
effects can be effectively described by a dynamical com- 
pact U(l) gauge field, and a SU(2) g aug e field that cou- 
ples to the color space of fi,A,a,a [28l. l29j. 

Instead of constructing the spin Hamiltonian, let us 
just consider a spin state, where the slave fermion fi t A,a,a 
fills a similar mean field band structure as the Fu-Kanc- 
Mele (FKM) model on the diamond lattice, and the flavor 
index a = 1, 2 plays the role of spin in the FKM model [5j : 

Ho = ~ tijfi,A,a,afj,A,a,a 

<i,j> ,A,a,a 

+ lX flA,a,Jab ■ (4 X 4)/j,A,fc,a. (18) 

The spin wave function is obtained after projecting the 
slave fermion mean field wave function to satisfy the 
gauge constraints: 

|G sp in) =I] P ^ = 2iV)®P(pf =0)|/ i>Aa , a ). (19) 



In Eq. I18[ when t%j is uniform and isotropic, there 
are three independent 3d Dirac fermions in the Brillouin 
zone; two of the three Dirac fermions can be trivially 
gapped out with an anisotropic , now we are left with 
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one Dirac fermion at Q = (0,0, ir), whose low energy 
Hamiltonian reads 

2N 

a=l A 

Ti — T z cr x , r 2 =rV, T 3 = T y . (20) 

Now we turn on a topological mass gap H4 to the mean 
field Dirac Hamiltonian: 

H 4 = ^2 m V'il,a r 4PAB^B,a) F 4 = T X . (21) 

With this mass, the band structure of the slave fermions 
with color index A — 1 is a topological insulator, thus 
A = 1 fermions have massless 2d Dirac fermion edge 
states; the slave fermion with color index A — 2 is a 
trivial band insulator without edge states. H4 also breaks 
the SU(2) gauge field down to another U(l) gauge field 
generated by p z , i.e. with nonzero H4, the slave fermion 
/ is coupled to two different compact U(l) gauge fields 
at low energy. 

Notice that in the current situation we did not partic- 
ularly assume the time-reversal symmetry in our model, 
thus one may expect the system to develop another mass 
gap H 5 — m^T^tp, where T 5 = t z o z . Indeed, in the 
original FKM model, this is precisely the mass gap that 
breaks the time-reversal and drives the system into a triv- 
ial insulator. However, the physical meaning of this mass 
gap is a slave fermion density modulation between two 
flavors and two sublattices, in the FKM model it is a 
spin density wave. Since our original spin model requires 
a constant slave fermion number on each flavor, then af- 
ter gauge projection this mass gap does not correspond 
to any physical order parameter in our lattice spin model. 

To build the connection with the NLSM Eq. HI we cou- 
ple the slave fermions with the antiferromagnetic Neel 
order parameter V: 



H 6 = m 6 iplX 5 ippP al 3 



(22) 



Both H4 and Hq are mass gaps of the Dirac fermion. 
Now the entire Lagrangian of the Dirac fermion can be 
written concisely as : 



U = cos((9)r 4 p z + sin(6»)r 5 -p. 



(23) 



IA is a unitary matrix. After integrating out the Dirac 
fermions, a 3+ld WZW term for unitary matrix U is 
generated [231 ] : 



WZW(W) = 



du 



d 3 xdr 



2tt 
480vr 3 



(U^dUf. (24) 



Once we assume there is a nonzero background H<±, the 
0— term in Eq. |H] can be derived by directly plugging 
U = cos((9)r 4 p z + sin(0)r 5 7> in this WZW term Eq. [Ml 



With a nonzero H4 in the bulk, the boundary of 
the system is described by 27V two-dimensional Dirac 
fermions with a SU(27V) symmetry, and these Dirac 
fermions are coupled to two U(l) gauge fields. The 
SU(iV) PCM at the boundary (Eq. [TO]) can also be di- 
rectly derived using the boundary Dirac fermions, once 
we break the SU(27V) symmetry to SU(A^)xSU(iV) x Z 2 . 
With the Z2 symmetry, the derived PCM model has pre- 
cisely 0' = 7T. 

If the Zi symmetry at the boundary is further bro- 
ken, a mass term can be turned on at the boundary: 
i?7 = myipflip. Just like the ordinary 3d topological in- 
sulator, this mass term 7777 drives the edge states into 
a quantum Hall state with Hall conductivity ±1/2 for 
the SU(iV)xSU(iV) charges respectively. With absence 
of topological degeneracy, a fractional Hall conductivity 
can only occur at the boundary of a 3d system. A do- 
main wall of 7777 is an analogue of domain wall of 0' in 
Fig. [It, and using the slave fermions it is straightforward 
to show that there are nonchiral gapless states localized 
at the domain wall of 777.7. 

So far we have ignored the dynamical gauge fields, 
which in 3+1 dimensional space-time can have a gap- 
less photon phase. In order to make sure the bulk is a 
fully gapped SPT, we need to drive the system into the 
confined phase of the U(l) gauge fields. Confinement of 
U(l) gauge field is driven by condensation of the mag- 
netic monopoles. In the free electron case, a monopole 
in a topological band insulator will carry gauge charge 
due to the topological 0— term in the electromagnetic re- 
sponse function [30j. The 9— term leads to "oblique con- 



finement" after the monopoles condense [3lJ, [32j. But 



in a system where charges are strongly interacting, the 
quantum number of the lightest monopole, as well as the 
nature of its confinement transition is not obvious, and 
I will leave this to future studies. Condensate of bound 
state between monopole and gauge charges in strongly 
interacting system is currently under study by other au- 
thors [33|, and here it is assumed that an ordinary con- 
fined phase is still possible by condensing appropriate 
bound state of monopole and gauge charges. The non- 
trivial edge physics of the 3d SPT will survive the confine- 
ment, because for example the 1+ld CFT at the domain 
wall in Fig.[TJ> cannot be gapped out without backscatter- 
ing between left and right moving modes, i.e. the domain 
wall CFT is always stable unless the SU(A r ) L xS\J(N) R 
symmetry is spontaneously broken down to its diagonal 
SU(A r ) subgroup. 

Just like the Id Haldane phase, the state described 
above is a 3d SPT only if the color singlet constraint 
fiP^fi — is strictly imposed on every site. By contrast, 
if this constraint is softened, namely the representation 
on every site is no longer the one in Fig. [TJ;, another 
mass term can be added to the mean field band struc- 
ture of the slave fermion: H$ = msip^sP 2 ^, which will 
completely destroy the bulk SPT, and gap out the edge 
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states without degeneracy. m§ corresponds to a "color 
density wave" on the lattice, which is not allowed with 
the on-site color singlet constraint. 

Our construction only applies to the self-conjugate rep- 
resentation in Fig.[T]:, which is invariant under the center 
Z 2 n of group SU(2iV), while the fundamental represen- 
tation is not invariant under the center Z 2 ^. Thus the 
3d SPT state constructed in this work cannot be classi- 
fied using the cohomology of SU(2iV) group, it might be 
classified using the cohomology of SU(2iV)/ Z 2 n- 



3.3 Edge states 

Let us take N = 2 as an example, and reinvestigate 
the domain wall configuration in Fig. [TJ>. Let us couple 
the SU(2)i charges to a U(l) gauge field A^a z , which 
is a spin gauge field that couples to a z only. Based on 
the gapless domain wall states, one can show that if a 
2ir— flux of is inserted at the origin x = 0, this flux 



11, 3 



111 . |12| . |34| , then inserting 



-flux of inside the domain wall will at least ac- 



will acquire gauge charge 2 

As a comparison, let us make a similar domain wall 
in a pure 2d system described by the same SU(2) PCM 
Eq. riOl with SU(2) L xSU(2) fl x Z 2 symmetry, and in- 
side the domain wall 9' = 27r, where the SU(2)l and 
SU(2)fl charges have opposite Hall conductivities. Since 
there is a Z 2 transformation connecting the two different 
sides of the domain wall, and the Z 2 symmetry exchanges 
SU(2)l and SU(2)h, thus the systems inside and outside 
the domain wall should have opposite Hall conductivities 
of A^ . Since for a 2d bosonic system without fractional- 
ization and topological degeneracy, the Hall conductivity 
can only be an even integer 
a 27T- 

cumulate 4 gauge charges. This proves that with the 
SU(2)ixSU(2)fl x Z 2 symmetry, the domain wall states 
at the edge of the 3d SPT described in Eq. |9] cannot be 
realized in a 2d system. 

Without the Z 2 symmetry that connects the two sides 
of the domain wall, the argument above would fail. 
For example, the 2+1 dimensional PCMs discussed in 
Ref. [3l have no such Z 2 symmetry that connects O' = 
2tt and 9' = 0. 

The same argument can be generalized to the case with 
N > 2. Now we conclude that Eq. M describes a 3d SPT 
protected (at least) by the subgroup SU(A f )xSU(7V)xZ 2 
of the SU(2iV) spin symmetry. Since SU(7V)xSU(iY)xZ 2 
is a subgroup of SU(2iV), the original SU(27V) invariant 
model Eq. [5] should also describe a 3d SPT. 



4. DISCUSSION AND SUMMARY 

In this work, we studied one class of 3d symmetry pro- 
tected topological phase, whose stability requires spin 
symmetry SU(2iV) or its subgroup SU(A^)xSU(A^)xZ 2 , 



but does not require other discrete symmetries. For 
large enough N, homotopy groups ^2d[ V (N)xv(N) ] and 
7i"2d-i [SU(iV)] are always Z for d > 2, thus our formalism 
and results can be generalized to any odd spatial dimen- 
sion. A 0— term can be defined for Neel order parame- 
ter V e u(jv)xu(jv) m an y spatial dimension. After 
breaking the SU(2iV) symmetry to SU(iV) xSU(A0 x Z 2 
symmetry, this bulk 0— term always reduces to a bound- 
ary e'-term with 6' = tt in the SU(iV) PCM at the 
boundary. 

This work is supported by the Alfred P. Sloan Foun- 
dation, Hellman Family Foundation, and NSF Grant No. 
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